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We study the quantum correction to conductivity on the surface of cubic topological Kondo 
insulators with multiple Dirac bands. We consider the model of time-reversal invariant disorder 
which induces the scattering of the electrons within the Dirac bands as well as between the bands. 
When only intraband scattering is present we find three long-range diffusion modes which lead 
to weak antilocalization correction to conductivity, which remains independent of the microscopic 
details such as Fermi velocities and relaxation times. Interband scattering gaps out two diffusion 
modes leaving only one long-range mode. We find that depending on the value of the phase coherence 
time, either three or only one long-range diffusion modes contribute to weak localization correction 
rendering the quantum correction to conductivity non-universal. We provide an interpretation for 
the results of the recent transport experiments on samarium hexaboride where weak antilocalization 
has been observed. 

PACS numbers: 72.15.Qm, 73.23.-b, 73.63.Kv, 75.20.Hr 


I. INTRODUCTION 

Samarium hexaboride (SmBg) along with PuBg and 
YbBg have recently emerged as prominent candidates 16 
for hosting topologically protected metallic surface 
states? 1 In particular, SrnBg - a material in which strong 
hybridization between samarium conduction d-electrons 
and strongly correlated /-electrons drives an onset of an 
insulating state at low temperature^^ _ has recently 
came into focus of theoretical and experimental studies 
as a most prominent candidate for the first correlated 
topological insulator. 13 " 23 

In order to experimentally establish the existence of 
the helical conduction time-reversal invariant states on 
a surface of a generic topological insulator using either 
transport or thermodynamic measurements, one needs to 
show that (i) the conduction at low temperatures is lim¬ 
ited to the surface; (ii) an inclusion of a small amount 
of magnetic impurities leads to localization of the sur¬ 
face states, i.e. presence of the time-reversal break¬ 
ing scattering potentials leads to localization; (iii) single 
particle states have linear momentum dispersion along 
the surface and no dispersion in the normal to the sur¬ 
face direction and (iv) there is a strong spin-orbit in¬ 
teraction, which leads to coupling between the momen¬ 
tum and spin of the conduction electron giving rise to 
the helicity of the carriers. In samarium hexaboride, 
a series of state-of-the-art transport studies have un¬ 
ambiguously shown that the resistivity plateau at tem¬ 
peratures be low 5 K elvirPEl is governed by surface con¬ 
duction onlyp3n2EH Furthermore, Kim, Xia and Fisk! 18 -! 
have examined the low-temperature transport proper¬ 
ties of the alloys Smi-xArBg for the non-magnetic yt¬ 


trium and ytterbium (A=Y,Yb) and magnetic gadolin¬ 
ium (A=Gd) substitutions. They found that while small 
amount (~ 3%) of gadolinium leads to insulating behav¬ 
ior in resistivity, substitutions of non-magnetic ions do 
not cause destruction of the metallic surface states.^ To 
verify that the conduction electrons on the surface have 
Dirac dispersion, G. Li et alphl have experimentally stud¬ 
ied the quantum oscillations of magnetization under ap¬ 
plied external magnetic field. By plotting the dependence 
of the index n which labelled the positions of the maxima 
in magnetization versus the inverse of magnetic held, G. 
Li et al. have shown that there is a contribution corre¬ 
sponding to the zero energy state which would only be 
possible for conducting state with Dirac-like dispersion. 

The magnitude of the spin-orbit interaction for the 
surface electrons in topological insulators can be indi¬ 
rectly probed by studying the quantum interference cor¬ 
rection to conductivity^ 428 upon decrease in tempera¬ 
ture, ST < 0, increase in conductivity (6a > 0) would 
signal weak anti-localization effect as opposed to weak 
localization corresponding to decrease in conductivity 
(<5cr < 0). The sign of the correction to conductivity is de¬ 
termined by the ratio of the spin-orbit scattering length, 
Iso, to the dephasing length 1$: for weak spin-orbit cou¬ 
pling, Iso l<f> (here l is a mean-free path) correction 
to conductivity is negative, while in the opposite limit of 
strong spin-orbit coupling Iso ^ l<t> and the interference 
correction to conductivity is positive. 

Thus, helicity of the Dirac-like carriers on the surface 
of topological insulators, associated with the fixed polar¬ 
ization of electron spin perpendicular to the momentum 
direction (i.e. Iso ~ l, here l is a mean-free path), neces¬ 
sarily leads to weak-antilocalization due additional Berry 
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FIG. 1: (Color online) Surface band structure for the strong 
topological insulator emerging from the inversion of the odd 
and even parity bands at the X points of the bulk Brillouin 
zone. Arrows denote the spin texture of the surface carriers 
corresponding to the groun d state configuration with the Tg 
quartet of the /-levels “12ii We assume that chemical potential 
crosses all three bands. Without loss of generality we consider 
the bands with the same chirality. In addition we will neglect 
the ellipticity of the Dirac pockets, but take into account the 
difference in the Fermi velocities of the electrons in different 
pockets. 


phase acquired by t he ca rriers as they scatter along the 
time-reversed paths. 25 29 

Application of an external magnetic field perpendic¬ 
ular to conducting surface destroys quantum interfer¬ 
ence processes leading to positive or negative magneto¬ 
conductivity - another signature of weak localization 
or weak anti-localization. The corresponding mag¬ 
netic field dependence of the conductivity correction is 
then described by famous Hikami-Larkin-Nagaoka (HLN) 
formula!^ 


Acthln(-B) = 


2tt 2 Tl 


l 0g * _ ^ ( I + ^ 
S B V l 2 B 


( 1 . 1 ) 


where a is a dimensionless parameter determined by the 
number of conduction channels and the strength of the 
spin-orbit coupling, ip is the digamma function, B 0 = 
H/4el Moreover, a > 0 for the case of the strong spin- 
orbit coupling and each independent conduction channel 
contributes 1/2 to the value of a, so that a = 1 for 
the case of Rashba-split bands on the surface. For the 
topological surface states, a = 1/2 for the case of a single 
Dirac band and a = 3/2 for the three Dirac bands. 

Recently, S. Thomas et alP^ have studied weak anti¬ 
localization effect and magneto-conductivity in samarium 
hexaboride. By fitt ing the experimental data with HLN 
formula, Eq. (1.1) for the case of a single band, the 


value of the parameter a came out to be approximately 



FIG. 2: (Color online) Schematic plot for the crossover behav¬ 
ior of the dimensionless coefficient a, appearing in Hikami- 
Larkin-Nagaoka e xpre ssion ( |1.1| ) , as a function of the ra¬ 
tio Bo/Bi Eq. (1.41 for the quantum correction to con¬ 
ductivity for the surface states in topological Kondo insu¬ 
lators. The value of the remaining parameter is chosen as 
Bo = 0.9Bi. Within the three-band model for the topological 
surface states, in the absence of the interband scattering there 
are three diffusion modes, which govern the quantum correc¬ 
tion to conductivity with a ~ 3/2. However, in the presence 
of the interband scattering processes two out of three diffusion 
modes become gapped, so that a ~ 1/2. However, at mod¬ 
erately high temperatures when becomes comparable to 
the gap of the remaining two diffusion modes also contribute 
leading to a ~ 3/2. Inset shows two independent scatter¬ 
ing processes from state with momentum k in band a to a 
state with momentum — k in band a' leading to the quantum 
correction to conductivity. 


equal to one, a ~ 1, for several sets of data. In the 
most recent transport experiments by Y. Nakajima et 
alPi also observe weak-antilocalization effect and their 
results seem to be generally in agreement with the earlier 
studies. On the other hand, these observations would 
contradict a natural expectation that the value of alpha 
must be close to a ss 3/2, whic h is due to the three Dirac 
bands on the surface of SmB a PEE2E31. 

Motivated by these experimental results, we calculate 
the quantum interference correction to conductivity in a 
generic cubic topological Kondo insulator. Our main goal 
is to account for the fairly wide distribution of values for 
the par amete r a obtained by analyzing the experimen¬ 
tal dataP^l As we argue in this paper, the presence 
of the disorder-induced interband scattering suppresses 
two antilocalization modes and reduces the value of a 
from naive a = [1/2 x number of zones]. We consider 
the surface band structure which consists of three Fermi 
pockets: one at surface T point and two at X and Y 
points of the two-dimensional Brillouin zone, Fig. [TJ We 
assume that the electron scattering within each band - 
intraband scattering - provides the strongest scattering 
mechanism, so that elastic scattering time is the shortest 
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time scale in the problem. Consequently, the disorder 
scattering between various bands is considered as a cor¬ 
rection to the intraband one. We show that (i) for the 
case when only intraband disorder time-reversal invariant 
scattering is present, all three conduction channels (per 
surface) will contribute to the interference correction to 
conductivity; (ii) the inclusion of the interband scattering 
shows that two conduction channels are suppressed while 
the remaining one contributes to weak anti-localization 
effect. However, at higher temperatures when the inverse 
dephasing time r^ 1 ~ T p (p > 0) becomes comparable 
with the size of the gap in the spectrum of the diffusion 
modes, we still find that all three modes contribute to 
conductivity. By considering as a parameter, we de¬ 
scribe this crossover behavior by showing the dependence 
of the parameter a on schematically on Fig. [2] For 
the analysis of magnetoconductivity using HLN formula 
(1.1) our result implies that the size of the correction for 


small fields and the large fields generally would corre¬ 
spond to diffe rent values of a. Specifically, we will show 
that the <0> generalizes to 


A " <B) = wn Y. 

i =0 


OLi 


log — - ^ ( 1 + — 
6 B M2 B 


, ( 1 - 2 ) 


where are the dimensionless parameters, which de¬ 
pend on the diffusion coefficients of the surface electrons 
from each band and scattering times and B\ 2 are de¬ 
termined by the interband scattering times, which give 
rise to the gap in the long-range diffusion modes, while 
Bi ~o = fi/Ael^ is still defined by the inelastic dephasing 
length Z^. In the limit of infinite interband scattering 
times - B 12 become equal to B 0 and we recover the HLN 
formula. In the limit of low magnetic field Bi/B 1 
(1.2) simplifies to 


A a(B ->• 0) « - 


247 l 2 Tl, 


x> 


i =0 


(1.3) 


If we now compare this expression with (1.1) in the limit 


of small magnetic fields, we see that parameter a can be 
written as 



We see that in the case when temperatures are not very 
low, Bq ~ Bi and even the gapped diffusion modes will 
nearly equally contribute to the localization correction. 
Therefore, depending on the temperature at which the 
experiments are performed and on the surface disorder 
one, expects that the values of the parameter a extracted 
by fitting the experimental data using may vary 

from a ~ 1/2 to a ~ 3/2 for the case of three Dirac 
bands. This situation is schematically shown on Fig. [2] 
Our paper is organized as follows. In Section II we 
introduce the model for the surface states. In Section III 
we present the calculation of the interference correction 


to conductivity for the uncorrelated mixture of the intra¬ 
band and interband disorder potentials. In Section IV 
we generalize the results from the previous Sections to 
the case of the correlated mixture of the scattering po¬ 
tentials. Sections V and VI are devoted to the discussion 
of our results and conclusions. Throughout the paper we 
adopt the energy units h = c = 1. 


II. SURFACE STATES IN THE PRESENCE OF 
DISORDER 

In this Section we setup the model and introduce the 
parameters which will be used in the calculation for the 
interference correction to conductivity. 


A. Hamiltonian and correlation functions 


The Hamiltonian for the surface electrons in cubic 
topological Kondo insulators can be written as a sum 
of three terms, which describe electron s near 1'. X and 
Y points in the 2D Brillouin zone, Fig. [ij 2 l 5 l 32 l 3 3 ] 

H= Y Y^jp^^'^pv C 2 - 1 ) 

j=F,X,Y per 


Here we neglect the anisotropies in velocitie s alo ng x and 
y-direction in X and Y pocketsP^Hln Eq. (2.1) momen¬ 
tum is taken relative to the center of the pocket. Intro¬ 
ducing the six component spinor 

'k i = [Vvpt i’Tpi VbfpT i’Xpl V'Vpt i’Ypl] (2-2) 


the Hamiltonian can be compactly written as follows 


/Cr 0 0 \ 

H = vY. ■ p, E = n„ <g> a, n„ = 0 Of 0 (2.3) 

\0 0 CyJ 

and the coefficients Cr,x,Y = Vt,x,y/v account for ve¬ 
locity anisotropies on different pockets. The underlying 
cubic symmetry requires Vx = vy, however there is no 
symmetry constrains on the ratio of velocities at T and 
X, Y pockets, so that generally r>r 7^ vx,y ■ In addition, 
we define the retarded and advanced Green’s functions: 


GZ’ A ( p,e) 
G$ A ( p,e) 


G?r A (p,e) 

G?/(p,e) 

0 

(e ± iS)a 0 + Vj(a ■ p) 

(e ± iS) 2 — v 2 p 2 


Goy A ( P.<0 


(2.4) 


Note that in the Hamiltonian (2.3) we ignore the higher 
order terms in momentum as well as other type of con¬ 
duction channels which may arise due to polarity driven 
bandspS strong surface potential 33 ^ etc. We will discuss 
how the presence of additional terms in the Hamiltonian 
may affect our results in Section [V] 
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B. Intraband disorder 


In the following we construct a theory of the metal¬ 
lic conductance of the disordered surface state of SmB 6 
using perturbative expansion in ppl 1. The values es¬ 
timated from experiments pfI ~ 100 suggest the pertur¬ 
bation theory to be a good approximation for the datcPR 
We consider TV, impurities with potential Uq and matrix 
structure described by U leading to the following expres¬ 
sion 

V^(r) = uqU^^ 5 (t - Rj). (2.5) 

i 

Averaging over an impurity ensemble yields 

b(p) = V«o^ p ' Ri , (V(p))dis = 0, 

( 2 . 6 ) 

(V r (pi)V r (p 2 ))di S = ruAulUU 5(p 1 + p 2 ), 

where A is the surface area and n, : = N^/A. In the 
following we first consider the intraband disorder, which 
implies that matrix U has the following block-diagonal 
structure: 


U = L 


3x3 1 


1 &0, 


(2.7) 


where / 3 X 3 is a unit matrix. We first evaluate the disor¬ 
der averaged correction to the self-energy: 

E(e) = n iU 2 J ^UG R ’ A (p, e)U (2.8) 


Clearly, E(e) has a block-diagonal matrix structure simi¬ 
lar to th e exp ression for the retarded and advanced prop¬ 
agators (2.4). It follows 


£j( e ) ~ Tj Q 1 = 7 vrnvjul. (2.9) 

2t j0 


where Tj o is elastic scattering rate, we neglected the real 
part since it leads to a small correction to e ~ VjPFj, PFj 
are the corresponding Fermi momenta and v ., is a single 
particle density of states per spin Vj = ppj/2-KVj. The 
corresponding expressions for the renormalized retarded 
and advanced correlators become 


Gf A 


(P, 


<0 


(e±i/2Tjo)ao+Vj(a ■ p) 

(e ± i/2T j0 ) 2 - vjp 2 


( 2 . 10 ) 


Similar expression for the correlation functions has been 
found for graphene.^ 


C. Classical conductivity 


Next, we can use the expressions for the correlators 
(2.10) to calculate the classical conductivity, which is 
given by 


c a p{u) = ^ 


J {v a G R (p,e + uj)v 0 G A (p,e)} . 


( 2 . 11 ) 


Here v a are the components of the velocity defined by 
the momentum derivative of the Hamiltonian (2.3) v a = 


V Pa H = uE„. The block-diagonal structure of matrices 


entering into (2.11) allows one to write 


t (0) 

a/3 


(w) « S a0 Y 


e v J v j T jO 

2(1 — icoTjo) 


( 2 . 12 ) 


where the subscript denotes that we have neglected the 
vertex corrections and on the last step we have assumed 
e {w, t~ 0 1 }. The vertex corrections can be formally 
included by making the following substitution in (2.11): 


A 


d 2 p 

w 


U6 A {e,p)v a 6 R (e + u,p)U (2.13) 


Summing the resulting geometric series to all orders we 
obtain 


e 2 VAV 2 T tI „ 

c a p{u) = Y 9n w TtI ' j = 2Tj0 ' ( 2 ‘ 14 ) 


Thus, we find for the Dirac electrons the transport life¬ 
time for the conducting states is twice the elastic scat¬ 
tering timeP^ 



FIG. 3: (Color online) Dyson equation for the Cooperon prop¬ 
agator V a / 3 n s(oj, q). Here the Greek indices encode both spin 
and Dirac cone - ’’valley” - components, solid lines represent 
the single particle propagators and U accounts for the disor¬ 
der potential. 


D. Cooperon Propagator 


The quantum interference correction to conductivity 
(2.14) is associated with the two-particle correlation func¬ 
tion known as the Cooperon. It satisfies the Bethe- 
Salpeter equation which is represented diagrammatically 
on Fig. [3] In what follows, we will use separate nota¬ 
tions for the ’’valley” and ’’spin” indices: we adopt latin 
superscripts for valley a = 1,2,3 and Greek subscripts 
for components of the Kramers doublet. Equation for 
the Cooperon propagator reads: 


r$S(".q) = Ao u%u% 

+ A o Y / 0 U i 6R } a a' C y(e + w, p) (2.15) 
X r$$( W) q )U%, [G%i(e, q - p), 
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with Ao = riiUQ and the summation goes over repeated 
indices. Since are diagonal in ’’valley” indices, we 
have 

[G^]^( e ,p) = [G^(e,p)] a ^ ab (2.16) 

Next for convenience we introduce the spin singlet and 
triplet components of the Cooperor l 35 * 36 ! and define 

q) = \ 

(2.17) 

where 5^2 = 0 ,x,y,z. Using the following identities 


( ^scr y ) a fi(cryas) flU — 28 a. (3 [ii 

S=0,x,y,z 

a=x,y,z 


(2.18) 


we can now express the components of Cooperon (2.15) 


in the right hand side in terms of matrices (2.17) using 


= )E (°S 1 <Ty)l3aCs b 1 's*{u,q L )(a y as 2 ) 7 6, 
Si,S 2 

(2.19) 

which follows directly from relations (2.171 and ( | 2 . 18 1 . 
To obtain the equation for the singlet and triplet com¬ 
ponents o f the Cooperon, we multiply both parts of the 
equation ( 2.15[ ) by the product of Pauli matrices intro¬ 
duced in (2.171). 


III. INTRABAND SCATTERING 


A. Gapless cooperon modes 


under the integral we recall the definition of the Greens 
functions ( 2 . 10 ) which are diagonal in the band index 


a, b. Electric c onduc tivity is given by a trace over band 
indexes in Eq. (2.11) and therefore in the absence of in¬ 


terband scattering only Cooperon components that are 
diagonal in band indexes contribute to quantum correc¬ 
tions to conductivity. For a diagonal Cooperon each term 
in Eq. (3.1) is diagonal in band indexes and therefore the 


system of equations splits into a set of independent equa¬ 
tions for each band. Thus, we only need to consider the 
components Cg iS 2 (w,q) defined in the same band, Fig. 
[4] After short calculation we find that the equation for 
the components of the Cooperon matrix in the valley 
a = r, A', Y can be compactly written as follows: 


M a ■ C a = A 0 / 4x4 . A q = r-Q 1 Ao 


(3.2) 


and the elements of the matrix M are given by Eq. (A9) 


in Appendix[A] In limit of u>r a o <C 1 and v a q <C 1 for the 
matrix we find: 


M n = 


-iu -o v aQx 


? VaQx 


(A-“) 


0 


2 V aQy 

0 

(A-“) 

0 


i 

TaO- 
(3 ' 3) 

To find the components of the Cooperon matrix we need 
to find an inverse of the matrix M a . The quantum correc¬ 
tion to conductivity is determined by the diagonal com¬ 
ponents C° A of the Cooperon matrix. The singlet compo¬ 
nent is given by 


CS 0 (u,q) = 


An 


^a T aO q 2 


■ *WT q0 


(3.4) 


while for the three triplet components we get 


To keep our calculations transparent, in this section 
we analyze the system with intraband scattering only 
and postpone treatment of the interband scattering for 
the next section. Here the matrix for the intraband dis¬ 
order potential is diagonal in both spin and valley in¬ 
dices fl2.7|) . In the momentum integral we insert equa¬ 
tion ( |2.19[ ). This gives an overall prefactor of (1/2 ) 2 in 
front of the second term, but the subsequent trace over 
the product of Pauli matrices will cancel one of them. It 
follows 

< -'S b 1 S 2 b (^iH) = A 0 8aa'Sbb'Ss 1 S2 
+ j ~^( a y rJ S 1 )a/3[Ga (e + w, p)]a 7 ' (3-1) 

X [Gjf (e, q — p)]pg'(as(Ty)s>jiCg b g^' b (w,q) 


C a xx {u,q) = A 0 - 


v W T l o(! + sin 2 </>,) - 2 iujT a0 


v " T aoq 2 


“a'aO'l * w7 a0 

C a = Au ^qg 2 ' r ao( 1 + C °S 2 <t>q) ~ 2 ^' r a0 

VV ’ ° ^a T a09 2 ^ *WT a0 

C* z {u,q) = Ao- 


(3.5) 


Thus we find that for the intraband disorder out of four 
Cooperon modes per each Fermi pocket, only one mode - 
singlet modes Cg 0 (uj, q) - remains gapless: propagation of 
electrons from Fermi pockets remains uncorrelated. This 
is not surprising given the strong spin-orbit coupling for 
the surface electrons. 


B. WAL correction to conductivity 


Clearly, there are two possibilities: (i) when the retarded 
and advanced propagators belong to the same valley, i.e. 
a = b = a' = b' and (ii) when they belong to different 
bands or valleys, a ^ a',b 7 ^ b'. To evaluate the trace 


As we have just demonstrated, the singlet channel is 
the most singular one (4.7), so that we can ignore the 
correction to conductivity arising from the triplet compo¬ 
nents. Thus, for the singular weak localization correction 
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(a) # 


for the magneto-conductivity 5<j(B): 



(b) • o 



^ = 2^ 



Bo 

B 


l(j>a 

XT 


- V]/ 


Bo 

B 


(3.8) 


where we took into account the leading (gapless) con¬ 
tribution to the Cooperon, l^ a = Z) a T 0 = l/4eB 0 , 
Z 2 = D a T a o and 'L(z) is the digamnra function, and 
we assume that is the same for all bands. The first 
digamma function in this expression can be replaced with 
log(Bol^ a /Bl^) due to the large value of its argument and 
the log (l^ a /la) drops out of Aa(B) = 5a(B) — <5ct(0). As 
a result, we can immediately identify the pre-factor a in 
Eq. ( |1.1| from Eq. (3.8) as a = 3/2. This value of a is 
a simple consequence of the fact that for a three band 
model without interband scattering each band gives a 
universal contribution to the WAL correction to the con¬ 
ductivity. 


IV. INTERBAND DISORDER 


FIG. 4: (Color online) (a) Pseudospin is coupled to momen¬ 
tum via strong spin-orbit coupling. As a result only one mode 
(I t-f) — | it})/\/2 corresponding to a total pseudospin 5 = 0 
state - singlet mode - survives; (b) Propagation of electrons 
from different bands is completely uncorrelated as bands have 
different Fermi energies and velocities (see text for details). 
As a result we only need to consider a = b case. Even in the 
case of the interband scattering, the propagation of the chiral 
electrons still remains coherent. 


(see Appendix B for details) we find 


Now we consider the disorder potential which also 
includes the component which induces the scattering 
between the different pockets: we expect the latter 
to remove the divergent nature of some of the singlet 
Cooperon components Coo- Without loss of generality 
we will consider perhaps the simplest type of the inter¬ 
band disorder potential: 

V(r) = Uo^U 0 S(f- fi) +u x y^U x S(f-fj), 

i i (4.1) 

Uo = To x (To, U x = T x X (To- 


(5(7(0) « e2p a,V 2 a T lo /^OO^.q) 

_r< v v J 


2 =r,x,y 


dq 

(2Xp 


(3.6) 


Carrying out the momentum integral, we find 


mo) = A E *(£). 

a=F,X,Y V aU/ 


(3.7) 


where is the phase coherence time. Since (5<r(0) > 0 
the inteference correction leads to weak antilocalization. 
Lastly, we remind the reader that on symmetry grounds 
Txo = t Y o Tro- 

In an external magnetic field, the momentum inte¬ 
gral should be replaced with the sum over Landau lev¬ 
els. Specifically, in the presence of the perpendicular 
magnetic field B , the momentum is quantized: q ^ = 
(n + l/2)Zg 2 with l\ = l/4eil. Setting — iuj = t^ 1 in the 
expression for the Cooperon and performing the summa¬ 
tion over Landau levels yields the following expression 


Here, T x is a 3 x 3 matrix whose diagonal elements all 
equal to zero, while off-diagonal elements equal to one: 
[T x \ab = (1 — S a b)- Just as in the case of the weak antilo¬ 
calization in graphenef 35: one can show that this disorder 
potential captures the essential physical features needed 
to describe the long-range diffusion modes. Moreover, we 
assume that any given impurity does not scatter electrons 
simultaneously between different pockets and within the 
same pocket - no correlations between the impurity scat¬ 
tering - so that for disorder average we have 

{V(fi)V(f 2 ))dis. = AoU 0 U 0 5(fi - r 2 ) 

„ „ (4.2) 

T A x U x U x 8(t\ ^ 2)1 

where Ao = riioUg and A^, = rii X u x with rii X being the 
concentration of the interband scatterers. We define the 
corresponding scattering times 

T ab = 7rn ixU x Vb, T a = T a0 + T ab , (4-3) 

b^a 
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with the second expression following directly from Eqs. 
( 2.8|4.1|4~2 |, while intraband scattering time r a o is de¬ 
fined by Eq. (2.91. The notations we adopt for the scat¬ 
tering times in (4.3) should be understood as follows: 
for a,b = T, X , Y The remaining scattering times are 
obtained by cycling permutation of these indices. Note, 
that our choice of the interband disoroder potential leads 
to Tpx = Tyx but at the same time Trx ^ Txr- The cor¬ 
responding self-energy correction to the Green’s function 
now reads 


2t„, 


(4.4) 


A. Cooperon modes in presence of interband 
scattering 

As we have already mentioned above our goal in this 
section is to verify if the gapless singlet Cooperon modes 
acquire a gap due to interband scattering effects. In what 
follows we will ignore the triplet Cooperon components - 
they are gapped already - and obtain the equation for the 
singlet Cooperon components only. The corresponding 
equation for the pseudospin components (in the valley 
space) of the Cooperon matrix can be obtained similarly 
to Eq. & It follows: 

C$’ a b (w,q) = Ao(5 aa '<5bf/ + A x [U x ] aa >[U x ]bb'+ 


bands a and 6 ; (iii) suppression of the interband interfer¬ 
ence Cooperon modes with a ^ b or c ^ d in Eq. (4.5) 
due to Fermi line assymetry. The latter effect is some¬ 
what analogous to the Fermi line trigonal warping effect 
in the band structure of graphene 35 ^ and is present even 
in the case of very small Fermi line missmatch with an 
important distinction that in SmB 6 considered here it 
suppresses interband interference. In SniBg it is likely 
that (p a — pb)l 1 for all bands and therefore the mech¬ 
anism (ii) is significant. Formally this means that in 
the case of interband terms in the Cooperon and Hikami 
boxes calculated in Appendix [A] and [B] would be smaller 
by factor of 1 /((p a — Pb)l as compared to the ones cor¬ 
responding to intraband terms. Therefore, we only need 
to consider nine elements in the Cooperon matrix: three 
diagonal ones Cq,Q , aa = Cq 0 and six off-diagonal ones: 
CQQ ,bb = Cqq. Consequently, we introduce the following 
notations: 


C. 


fcEo 

r-irx 

L/ oo 

°00 


r x 

L/ oo 

r<XY 

°00 

Wo T 

r'YX 

°00 

r<Y 

°00 


(4.7) 


a. Eigenvalues of the Cooperon matrix. To solve the 
equation (|4.5| ) we first consider the following eigenvalue 
problem : 37 


= 


E 

a'b' 


M g 


, 6 V ( 0 , 0 )^ 


(4.8) 


+ • A ^aa",fch"(w,q)C'o 0 6 b (w,q), 


(4.5) where the components of the matrix jCi are given by 


a"b" 


shows that M. a 


An 


A 4 aa ' bb' fa: q) n ^ ] 


( |4.6[ ) . The quick calculation (see Appendix A for details) 
66 '( 0 , 0 ) = Maa’SabSa'b ' with 


(4.9) 




/ JT_ 

rx 


r d 2 p 

M = | 

f r ro 
r r 

r rx 
rx 

rrv \ 

T Y I 


1 r xr 
l _ZE_ 

rx 0 
rx 

rxv 1 

TY_ I 


\ty r 

r yx 

r Y o / 


x [Gf(e + w,p)] Q: 7 /[Gb l (e,q- p)]ps>(vy)6'y8 a a'5bb' + 

+ E j {Vy)ap[Ga'(t + W, P)]cry' 

X fa*b' 91 P)]/56' fay')S' 7 ' fa,]bb' 

(4.6) 

for convenience and summations are performed over re¬ 
peated spin and pseudospin indices. 

In presence of inter-band scattering it is important to 
take into account the off-diagonal elements of the Cooper- 
ons in the band space, i.e. consider the full 9x9 matrix 
(with spin indexes S = Si = 0). This consideration 
is significantly simplified due to the effect of the Fermi 
line missmatch between different bands, due to T and 
X, Y bands having different Fermi velocities. Further¬ 
more, X and Y bands are characterized by asymmetric 
Fermi lines, see Fig. [T] The resulting missmatch of the 
phases of wave functions results in: (i) suppression of the 
disorder induced interband scattering matrix element due 
to the Fermi wavelength missmatch; (ii) suppression of 
the contribution of the interband terms to the conductiv¬ 
ity oc 1/( fa a —Pb)l), where p a ,Pb are Fermi momenta for 


The eigenvalues (4.8) of this matrix can be most com¬ 


pactly written in terms of the following parameters 


„ _ 1 / T a n \ _ 1 / T a 

r ab — 0 ( ta^b ) ■> ta — 0 I ^ 

^ \T~ab T~ba ) ^ \^~a0 


7± 


= - E ta ± 


\ 


E*« 


E r “ b 

a^b 


(4.10) 


with t a = 57 ^ — Note that since the parameters r a b < 
0 and t a <0 the parameters 7 7 ar e positive. Then the 
expressions for the eigenvalues (4.8) are 


Aq — 1, Xi : 2 — 1 — 7±> 


(4.11) 


Clearly, if we neglect the interband scattering, we find 
a threefold degenerate eigenvalue A = 1. As we have 
discussed above this situation corresponds to the exis¬ 
tence of three gapless modes. Inclusion of the interband 
scattering processes lifts the degeneracy leaving only one 
gapless mode. Below, we first show that j± determine 
the gap for the diffusion modes and then compute the 
diffusion coefficient for the gapless mode. 















b. Eigenvectors. From the analysis of the equation 
(4.16) it is clear that the diverging contribution to the 


Cooperon emerges for the eigenvalue Ao = 1 with the 
corresponding components of an eigenvector 


*$ = No S -?, 


(4.12) 


where a = T, X 1 Y and proportionality coefficient N 0 is 
the normalization constant. Similarly, the eigenvector 
components for an eigenvalue A = 1 — 7 + are 


= 


JVi 


(1 + 7 a )r 0 


- E 


NqNi 


(1 + 7 d)rl 




(o) 

ab ' 


(4.13) 


Here N\ is a normalization constant and 7 a = 
7 +t“ 1 /(t i 7£ + TjTy + t Yp). Similar expression can be 
obtained for the third eigenvector ^t ab '. We note, that 
the eigenvectors satisfy the orthonormalization con¬ 
dition, 


/ j ^ ab ^ ab ' 

a,b 


(4.14) 


Having computed the eigenvect ors, w e can now express 
the components of the Cooperon (4.8): 


c oo’ c 'V> q) = E Ai ( w > q)^ib^id> ( 4 - 15 ) 


Using this equation together with the normalization con¬ 
dition (4.14|, we can re-write the equation for the com¬ 


ponents of the Cooperon matrix (4.5) as follows: 

^ [(1 - A i)6 ik + W lfc (w, q)] A hj (w, q) = V iJt (4.16) 

k 

where the matrix elements W,:fe(w, q) are given by 
W ifc (w,q) = 

= ^ [M aa > l bb , (0, 0) — Maa',bb'(w, q)]'F^,\I/ 


(fe) 

b 


ab,a'b' 


(4.17) 


and they are obviously vanishing for q —>■ 0 and w —> 0 . 

The matrix elements V, 7 are obtained from expanding 

the first two terms in the right hand side of the Eq. (4.5) 

(i) - 

in terms of the eigenvectors 4 , ' ) / : 

Ao SacSbd + A x [U x ]ac[Ux\bd = E V ^ab^cl ( 4 - 18 ) 


It is straightforward to compute the coefficien ts Ag (to, q) 
by solving the system of linear equations (4.16). For 


example, the coefficients Aoj(cv. q), which contribute to 
the gapless Cooperon mode can be found by using the 
fact that the matrix elements Wik(u>, q) are small for 
small oj and q. It follows 


A 0j {ijJ, q_) 


V, 


oj 


W 00 (w,q)' 


(4.19) 


The matrix elements entering into this expression can 
be found from ( 4.18|4.17 ) using the normalization con¬ 
dition (14. 141). However, for our subsequent analysis of 


the quantum correction to conductivity we will also need 
the remaining two contributions to the Cooperon. After 
some algebra we find that the resulting expression for the 
diagonal components of the Cooperon ( 4.7|4.15 ) can be 
written as follows: 


cs « = 


Win 


where 


r “ ^ d q 2 - iu + IV 


1 = E 1 

n T a 


(4.20) 


(4.21) 


and the gaps in the denominator are 

r 0a = o, ri a = r a _1 7 + , r 2a = r ( 7 1 7_ (4.22) 

and Wi ~ 0 ( 1 ) are dimensionless parameters determined 
by the combination of the intra- and interband scat¬ 
tering rates. For reader’s convenience, on Fig. [5] we 
plot the dependence of the coefficients Wir and u; 2 p and 
the eigenvalues A.; as the function of the ratio rp/rro 
for the specific choice of the scattering times such that 

'Ca/'^ba “ 'Cb/'Cab- 



FIG. 5: (Color online) Plot of the eigenvalues A; (i = 0,1,2) 
and the weight coefficients which appear in Eq. (4.201 as a 


function of rr/rro- Here without loss of generality we con¬ 
sider a special case when the following relation between the 
scattering times holds: (a ^ b) 


In the expression for the singlet Cooperon we intro¬ 
duced the diffusion coefficient for the singlet long-range 
mode: 


D = H 


n 

T~aO 


D n 


A, 

Ta “ T ab 
b^=a 


(4.23) 
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Thus, we find that pseudospin symmetry breaking (band 
mixing) perturbations produce the relaxation in the 
spin-singlet components of the Cooperon except for the 
single mode, which is protected by the time-reversal 
symmetry.^ The trajectories in Fig. giving rise to 
the gapless Cooperon mode correspond to each ballistic 
segment the blue and red lines (moving in opposite direc¬ 
tions shown by arrows) reside in the same electron band 
in the BZ and combine to form a spin-singlet. This is 
the pair of trajectories the interference between which 
is protected by time reversal, since time reversal maps 
each band on itself. Note that this is in contrast to the 
case of graphene where the time-reversal symmetry maps 
the two different bands on each other and therefore the 
band-singlet Cooperon mode is protected by time rever¬ 
sal symmetry^. In the following Section we evaluate the 
quantum correction to conductivity appearing due to the 
presence of this single long-range diffusion mode focusing 
specifically on the value of the pre-factor a appearing in 
the HLN formula (1.1). 


B. WAL correction to conductivity 


conductivity 6a(B): 


Sa(B) = 


2i r 2 


E 


WOaD a T t 


4- 


Bo 

B 


2DT a0 

2 " 


2 - — ) 
TaoJ 


- 4- 


Bo 

B 


(4.26) 


where we took into account the leading (gapless) con¬ 
tribution to the Cooperon, = Dt$ = 1/4 eB 0 , 1% = 
Dr t and H^z) is the digamma function. The first 
digamma function in this expression can be replaced with 
log(.Bo R/Bl%) for l^/lt 1 and we recover the structure 


of the HLN expression (1.1) for Aa(B) = 5cr(B) — 5a( 0). 


Therefore, from our result for the correction to conduc¬ 


tivity (4.26) we can immediately identify the pre-factor 


in front of the square bracket with the the weight factor 
a in Eq. (1.1): 


Wpa-Dq 

2D 


TaOj T a0 


(4.27) 


We again disregard the correction to conductivity aris¬ 
ing from the triplet components as they are suppressed by 
the intraband scattering with the largest gap scale in our 
model. Moreover, we can omit gapped singlet modes due 
to the interband scattering keeping only the gapless mode 
identified in the previous section. This gapless mode is 
protected from dephasing due to the disorder scattering, 
but is still suppressed on length scales = yjD t$ due 
to the inelastic electron scattering characterized by de¬ 
coherence time Tt/,. We have, see Appendix B for details: 


For the moderately large magnetic fields, however, the 
remaining two modes will contribute to conductivity. 
Their c ontrib ution is formally given by the same expres¬ 
sion as (4.26) where we have to replace t^ 1 with r )±r~ l . 
It follows 


A<j gap (B) 



i —1 


OWia 

WOa 


X 





(4.28) 


5a ~ E 2 —-) A A 3 x 

a=r,x,Y' Tq0 ' 


x J Cqq(lo, q) - 


dq 


(4.24) 


( 2 tt) 2 - 

Carrying out the momentum integral, we find 


<5ct(0) 




2 

X ^2 Wia lo § 
i =0 


( n' 

Y max {r^ , 



(4.25) 


where r $ is the phase coherence time and we assumed 
r f _1 r )±r~ 1 (i = 1,2). Since <5 ct( 0) > 0 the inteference 
correction leads to WAL. 

Similarly to the calculation in Section [IllB] in presence 
of an external magnetic field the summation over Landau 
levels yields the following expression for the magneto- 


with 4eI3j = ( 7 ±/D) • max{r r 1 , t x 1 , t^ 1 }, so that the 


r- 1 r -i t -A 
r r > T x > t y 

total correction to conduct ivity becomes Adtot = Act + 
Act, 


gap- 


Note, that unlike (4.26), the contribution from 
AfXgap is temperature independent. 

As we have already discussed above, in the ab senc e of 
the interband scattering, the eigenvalue problem (4.8) be¬ 
comes degenerate, since in that case the second term on 
the right hand side of that equation vanishes, while the 
integral in the first term equals one. Therefore, in that 
case there will be three independent singlet long-range 
modes and in that case a = 3/2. For finite albeit small 
interband scattering we find a ~ 1/2 from Eq. (4.27). We 
note that the value of a is not universal due to asymme¬ 
try between of parameters for T and X(Y) pockets. This 
non-universality resembles non-universal weak localiza¬ 
tion correction in a metal with partially polarized mag¬ 
netic impurities when the rotational symmetry is bro¬ 
ken.^ 8 It is important, however, to keep in mind that as 
temperature is decreased, one may expect a crossover be¬ 


havior from a ~ 3/2, when 


4 eDB^ to a 


1/2 for 


-1 


4 eDB 7 , see Fig. 2. 
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FIG. 6: (Color online) Plot of the total correction to conduc¬ 
tivity A(j to t) as a function of magnetic field for a = 0.45. 
The values of the scattering rates and diffusion constants 
have been chosen as follows: Try = 1.2rrx, txy = tty and 
D x = D y = 0.95 D r . 


C. WAL in correlated disorder 


It will be instructive to consider the model for the mix¬ 
ture of correlated disorder, when the disorder potential 
at a given impurity site can scatter within the same val¬ 
ley as well as between different valleys, and compare the 
corresponding value of the weight factor in this case with 
the one found above (4.271. 


For the correlated disorder mixture we consider the 
following model: 


V(f) = J2 uU5{r- Ri), U = U 0 + (U x 


(4-: 


where Uq, U x are defined in Eq. (4.1) and coefficient 
C = u x /u. Now in addition to the two scattering times 
r a0 and r a defined previously ( 2.9|4.3 | we also introduce 
the relaxation time due to disorder correlations: 


The analysis of the Cooperon eigenvalues can now be 
done along the same lines as above. Specifically, the right 
hand side of the equation (4.8) will now acquire an extra 
term proportional to t~ x . Setting the eigenvalue Ao = 
1 and solving for the eigenvectors within the required 
accuracy yields: 




(o) 

ab 


= ^<U + (i-<U)( — ^ + — r — 

Ta. \ T ax T a Tjj X 7*6 


(4.32) 


where f* ~ r* appears as a result of normalization of 
the eigenvectors. An explicit expression for it in terms 
of scattering times is quite cumbersome and will not be 
given here. In addition, we also introduced 


i 31 

II 

1 - 

7”inter 

1 to 

k 3 

T a 0 

Tax 

T- 1 ~ 

inter 

T rx 

+ t xy 

+ Ty 


b^a 


(4.33) 


and rT t ' er = + 7y/ stands for the average 

interband scattering rate. 

From our results for the eigenvectors it is clear that the 
finite r ax scattering time introduces off-diagonal in val¬ 
ley indices components of the Cooperon matrix, which 
give non-zero contribution to conductivity. However, the 
contribution of these off-diagonal terms is sub-leading in 
power of T a jT ax to the diagonal ones. If the off-diagonal 
Cooperon elements are neglected, the resulting expres¬ 
sion for the coefficient a in the HLN formula reads 


= V^ 2-— —, (4.34) 

y 2D V W T a0 


where = v„t, 


and D is found by the same expression 
as in (4.23) by replacing Tq a with 7 \ a and r t -1 = fjT 1 + 

further reduce 


+ 7% 


=T 


r * 1 

the weak-antilocalization effect. 


Y 


As we expected, the finite r a 


V. DISCUSSION AND CONCLUSIONS 


Tax = TniC,u 2 v a = CtJ. (4.30) 

In turn, the self-energy matrix becomes off-diagonal in 
band indices due to disorder correlations: 


^R,A 
^ab 


(e) « T 


ie o 

2 r a 


T 


iU x 

2 7" 

^' an 


(4.31) 


As a consequence of this, the matrix Green’s function 
G -1 = Gq 1 — E also becomes non-diagonal in band in¬ 
dices. However, in the computation of the Cooperon ma¬ 
trix within the required accuracy we can neglect the pres¬ 
ence of the t~ x in the Green’s functions. The reason is 
that in the calculation of the singlet components of the 
Cooperon, the matrix elements which are proportional to 
t~ x are of the order of (pfO” 1 ^ 1 an d therefore can be 
neglected. 


Perhaps one of the main challenges in identifying the 
nature of the conducting surface states, such as the 
dispersion and helicity, in topological Hondo insulator 
such as SmBg lies in establishing to what extent the 
Dirac surface states remain well defined despite the fact 
that (a) hybridization between the conduction d- and 
/-orbitals is significantly reduced on the surfac^H an d 
(b) the band bending effects due to disorder scatter¬ 
ing on the surface, which leads to an appearance of 
the states inside the hybridization gap. 33 -l For example, 
recent ARPES measurements (see Ref. [40] and refer¬ 
ences therein) and low-frequency, radio-frequency and 
microwave conductivity® 3 seem to be in support of the 
picture in which conventional Rashba-split bands dom¬ 
inate the low-temperature transport properties on the 
surface of SmBg. On the other hand, recent radiation 
spectroscopy measurements and magneto-thermoelectric 
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transport results and well as spin-resolved ARPElPH sup¬ 
port the picture of the topologically protected surface 
states. Interestingly, the Nernst effect data on the (Oil) 
plane reports the effective mass for the carries of the order 
of 100 of bare electron mass in agreement with existing 
theoretical estimates! 32 * 33 * 

Observation of the weak anti-localization correction to 
conductivity in topological insulators generally serves as 
an indication of the strong-spin orbit coupling and, there¬ 
fore, is used to confirm the helicity of the conducting 
surface states. In topological Kondo insulators, how¬ 
ever, the analysis of the experimental data is complicated 
by the possibility of the conventional polar bands which 
will be split by the Bychkov-Rashba spin-orbit coupling 
A so- When disorder does not induce scattering between 
the Dirac and parabolic bands, one may expect a weak- 
antilocalization correction to conductivity provided the 
spin-orbit coupling is strong enough: 37 -1 Furthermore, 
the correction appears to be non-universal and is propor¬ 
tional to (\soPFT~tr)- 2 where r tr is a transport timeP^l 
Clearly, whether the magnitude of this correction is of the 
same order as the one we find for Dirac electrons depends 
on the magnitude of the A so'- at strong Bychkov-Rashba 
splitting of the parabolic bands their contribution to con¬ 
ductivity quantum correction may be strongly suppressed 
by this additional factor. Lastly, for the scattering which 
mixes the Dirac bands with the parabolic bands, one sin¬ 
gle diffusion mode is expected to be present leading to 
a ~ 1/2 for a given surface. 

Another correction to our results above may appear 
due to the presence of magnetic scattering on the sur¬ 
face, which may change the picture of weak localization 
corrections presented above. In fact, based on the mag¬ 
netoresistance data^S, it has been recently argued that 
unscreened /-electrons give rise to ferromagnetic state 
on the surface of SmBg. One expects therefore, that 
the spin-flip scattering on the surface will gap out all 
Dirac bands completely suppressing transport. Even in 
the case when only one Dirac pocket is not gapped, the 


quantum interference correction to conductivity will be 
strongly suppressed resulting in the values of a < 1/2. 
In addition, as it has been shown recently, the opening of 
the gap in the Dirac bands may actually change the sign 
in the quantum correction to conductivity from weak- 
antilocalization to weak localization. 2| 3 Perhaps the fact 
that this crossover has not been observed-D suggests that 
at least one of the Dirac bands remains ungapped. 

In this paper, we have considered the quantum cor¬ 
rection to conductivity in the model with three Dirac 
bands with intraband and interband disorder. We find 
that for the case of the interband disorder there is only 
one singlet long-range mode leading to the quantum cor¬ 
rection to conductivity. The resulting expression for the 
weight factor shows that depending on the ratio between 
the diffusion coefficients and interband scattering times, 
one may expect the smooth crossover from a ~ 3/2 to 
a ~ 1/2. In contrast with the single band case when 
a = 1/2, the presence of the interband scattering reduces 
the value of a so that it becomes non-universal. In fact, 
the interband scattering itself may be asymmetric, due 
to the large Fermi line missmatch between T and X, Y 
bands the scattering between them could be suppressed. 
In this case, independent contribution of the F in addi¬ 
tion to the mixed X and Y components would result in 
a ss 1 which may explain the value observed in Ref. l30l 
Our results are generally in agree ment with recent low- 
temperature transport experiment/ 30 31 on SmBg. 
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Appendix A: Components of the matrix M 


In this Section we provide an explicit derivation for the elements of the matrix M which enters into the equation 
for the Cooperon (3.21. This matrix is defined as follows 


1 

2 



{[Ga (e + w, p)] T (a y a Sl ) [Gf{e, q - p)] (cr s a y ) j 


M&siv, q) 


(Al) 
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First, we compute the diagonal components in pseudospin space: 

(e + u + i/2T a )(e - i/2r b ) + v a v b p 2 - v a v b pq 


Mf 0 >,q) = / 


[(e + w + */2r a ) 2 - v 2 p 2 ][(e - i/2r b ) 2 -u 2 (q-p) 2 ] J ’ 


q) = 


d 2 p / (e + w + i/2r a )(e - V 2r &) + [p w ~ Py) ~ Px(Qx ~ Px)] 


= I ^ 




4tt 2 ( [(e + w + i/2r a ) 2 — v 2 p 2 ][(e — i/2T b ) 2 — u 2 (q — p) 2 ] 

d 2 p f (e + a; + i/2T a )(e - i/2r b ) + -p x ) - p y {g y - Py)} 

4tt 2 \ [(e + w + i/2r a ) 2 — f 2 p 2 ] [(e — i/2r b ) 2 — r> 2 (q — p) 2 ] 

2 ~ ( (e + oj + i/2T a )(e~i/2T b ) -v a v b p 2 +v a v b pq 


(A2) 


[(e + w + i/2r a ) 2 - c 2 p 2 ][(e - */ 2r b) 2 - y l (q - p) 2 ] J ' 

We calculate each of these integrals separately. In what follows I use the following approximation 


wT a , b < 1, Vj\q-p\xiVjp 1- 


<7 cos 0 q 2 sin 2 0 


PFj 


+ 


2p 2 


Fj 


Consider the following integral 

dab = 


d 2 p 


1 


47T 2 \[(e + uj + i/2r a ) 2 - vlp 2 ][(e - i/2r b ) 2 - v%p 2 } 
and introduce the following variables £ = sign(e)v a p — e, r ab = v a /v b . We have 

sign(e)oo 

1 

lab = 


(A3) 


(A4) 


2 FVI 


[sign(e)^ + |e|]sign(e)ci£ 


[(e + uj + i/2r a ) 2 - v 2 p 2 ][rl b (e - i/2r b ) 2 - v 2 p 2 ] 

+oo 


dd 


(A5) 


2irvl 2e 2 (l + r ab ) J (w + */2r a - £)[( r ab - l)e - i/ 2r b ~ £] 


27TU b T ab 


1 - 2iuT ab - 2*(1 - r ab )eT ab ’ 


T ab= T a * l +T b X . 


Clearly, from this expression it follows that all off-diagonal elements of the Cooperon matrix will remain finite in zero 
momentum and frequency limit. Thus we need to analyze (A2) for a = b only. We have 


Mi 


00 - 1 + lUJTa - 


2 2 2 
y a<l T a 




cu+i 

- < 1 , 




1 + ibJTa — 


2 2 v a p Fa 

v 2 a q 2 rl v 2 q 2 r 2 cos(2ip q ) 




1 + ibJTa - 


2 2 2 


V a < f T a C0S(2 ip q ) 


2 4 

where we used 

TTU a T a A = 1. 

Next, we consider the off-diagonal components is pseudospin space. We have: 

d 2 p f (e + w + i/2r a )(q x - p x ) + (e - i/2r b )p, 


(A6) 


(A7) 


M ab (m a) = 

01[ ,qj J 4* 2 \[{e + u + i/2r a ) 2 -vlp?}[{e-i/2T b ) 2 -v 2 {q-p) 2 } 
M o|( w q) = j^P j (e + w + i/2T a )(q y - p y ) + (e - 


^12 (w,*0 = - 


[(e + w + i/2r a ) 2 - v 2 p 2 ][(e - i/2r b ) 2 - t; 2 (q - p) 2 ] 
d 2 p f (e + uj + i/2T a )p y (q x - p x ) + (e - i/2r b )(q y - p y )p x 


(AS) 


4tt 2 ( [(e + w + i/2r a ) 2 - , c 2 p 2 ][(e - */ 2r &) 2 ~ ^(q- P) 2 ] 


while the remaining components will give zero. Thus, collecting all the terms we obtain: 


M = 


2 ° - lu}T a ~% V aVTa COS <p q 

( 2 2 2 " " 

1 • . V n q r n . 

1 - IbJTa + °^ 2 + 4 

-5 v a9 2 T- 2 sin(2^) 


■fVo^T-oSm 

0 


0 


~%v a qT a sm p q 
l y lQ 2T isin(2p q ) 

( 222 222 /o \ 

1 - ibJTa + - Vaq T °4° ’ ) 0 

0 


(A9) 
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This result shows that only single mode corresponding to the singlet component of the Cooperon matrix will remain 
gapless. 


Appendix B: Quantum corrections to conductivity 


1. bare Hikami box 




FIG. 7: (Color online) Panels (a)-(d): diagrams contributing to the weak localization correction to conductivity. Panels (e)-(h): 
the same contributions as (a)-(d), shown in the representation of Hikami boxes. 


Contribution to conductivity from the bare Hikami box, Fig. [7](a,e), is: 


x (b e 
= 2 tt 


/ (k, e + w)] Q/j r^(w,q)[Gf (q - k, e + (q - k, e)| 7 i (Bl) 


To evaluate this correction we employ Eq. (2.191. Calculation of the trace over the pseudospin degrees of freedom is 
done by Mathematica. The resulting expression can be simplified by neglecting the dependence on external momentum 
q in the single particle correlators. In addition, as it follows from the calculation of the traces we can also neglect the 
frequency dependence in the denominators. We are interested in find the contribution from the most singular terms 
in the Cooperon. For the diagonal components of conductivity it follows 


E 


2 2 
e vt 


i=r,x,y 


dkdq 

J 


E 9 C s 1 s 2 (q) Tr {&a ( k , z)cnG*{ k, e + ujcrfa^ x 

x [Ga (q - k, z)\ T °i[Ga (q - k, e + w)] T <r y as 2 


SiS 2 


} 


To calculate the trace we use the following relations: 


r/V-R An \iT &y\tR,A& 0 {kx@x T kyCTy )]^"y T rn c 11 /■ n \ 

[G"’ A (k,e)] J = -j- 212 a i = [ 2 ^,0 - 1 ]<Ty<Ti<Ty, (t = 0 ,X, V ,z). 

e R,A V a K 

Using these relations and neglecting the q dependence in the nominators of the G A ' R for the trace we find 


1 [dt 

2 S Tr 


{G A (k, e)<JiG*( k, e + ur)af cr^ [G A (q - k, e)] T erf [G* (q - k, e + w)] T a y a S 2 
^SiS 2 e4 {^Si.o - 3(5 Slji - <5 Sl yj 


} 


where on the last step I used w <C e and I have also assumed v a k ~ e. Thus (B21 becomes 

v 2 a f dq 


E 


a=r,X,Y 


27r J (27t) 2 


[4G 0 a 0 (q) - 3G“(q) - Gff(q)] 


kdk 


2 tt [(e R + w) 2 - v 2 k 2 ) 2 [e 2 A - v%k 2 ] 


(B2) 


(B3) 


(B4) 


(B5) 
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We deal with the momentum integral as follows: 

oo sign(e)oo 

J f ( vk ) = ^~2 J {sign(e)e fc + |e|} sign(e)cfe fc / ((e + £ fe )sign(e)) 


— e 
+oo 


27 TV 2 


+ OO +OO 

J f(Z+\e\)dt = v J /(<£+|e|)d£ 


Finally, the result is 


E 


e 2 VaV 2 a Tl 


a=r,X,Y 


dq 

8(1 — iwr) 3 J (27r) 5 


[4C 0 “ 0 (q) - 3C“(q) - C£(q)] . 


(B6) 


(B7) 


2. first disorder correction to the Hikami box 

The expression for the second correction to conductivity, shown on Fig. B b ,f), reads 

5 <r\? =A 0 E ^ CI SiS2(q) Tr {^a (k,e)o- i Gf(k,e + cu)17G , f(p,e + a;)gJg^ 1 x 


i=r,X,Y 
VA 


(B8) 


X [Ga (q - P, e)] T o-f [Ga (q - P, e + u)] T U T [G* (q - k, e + w)] T cr y CTs 2 } . 


Here we took into account that only diagonal part of the disorder potential contributes to the conductivity correction, 
since correlation functions are diagonal in valley indices. Taking into account the expressions for the propagators and 
relations ( |B3[ ) we obtain 

e 2 Vg f dq 
2lT / 47T 2 


^ 2 ) =Ao E Et 5E KU(q) 


a=r,x,Y 
’ pdp 


SiS 2 


kdk 


J 2 tt [( £jj + w + ±) 2 - v 2 a p 2 } 2 [{e - £) 2 - 77 2 p2] 7 2 tt [(e fl + w + ^) 2 - u2/c 2 ] 2 [(e - £) 2 - v 2 a k 2 } 
x TV {[e A <Jo + u a ( k ’ S)\(Ji[{e R + w)a 0 + u a (k ■ <r)][(ej? + u)a 0 + v a (p • a)\cr Sl x 
x [cacto + u Q (p ■ a)\<Ji[(e R + w)<r 0 + ^a(p • v)][(zr + w)cr 0 + w(k • n)]crs 2 } . 

Computation of the trace and subsequent integration over momenta in the limit w = 0 yields 


(B9) 


a=T,X,Y 


e 2 v a vlT% ( r a 


dq 


[Coo(q) - c“(q)] ■ 


,T a0 J J (27r) 2 

Note, that an additional pre-factor appears since the relation time differs from the intra-pocket scattering time 


(BIO) 


3. second disorder correction to the Hikami box 


The third correction, Fig. [7jc,g), is the same as the first correction (BIO) to the Hikami box diagram: 

(3) _ e2v aV 2 a Tl f r a \ f dq 


= - E 

a=T,X,Y 




(27r) 


2 [ c oo(q) - c'S(q)] ■ 


(Bll) 


Finally, the fourth and the last correction to conductivity, Fig. [7](d,g), is small in para meter 1/pfI 1 and can be 
ignored. Adding up all three contributions to the conductivity we find expression (4.24) in the main text. 
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